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^NJ ! Abstract 
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^sO ' We propose the deterministic rate equation of three-species in the 

reaction - diffusion system. For this case, our purpose is to carry out 
\ the decay process in our three-species reaction-diffusion model of the 

form A + B + C — > D. The the particle density and the global reaction 
rate are also shown analytically and numerically on a two - dimensional 
square lattice with the periodic boundary conditions. Especially, the 
crossover of the global reaction rate is discussed in both early-time 
and long-time regimes. 
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Recently, there has been considerable interest in the front propa- 
gation under the initial segregation of two - species reactants in early 

- time reaction - diffusion process. It is furthermore well - known 
that the binary reaction has mainly been investigated on the process 
of ternary reactions 1 ) and many phenomena in nature. 2_6 - ) The early 

- time behavior for reaction - diffusion system have been argued by 
Taitelbaum et al, 7 ) and it has precisely been shown that both the 
global reaction rate and the reaction front increase as a function of 
t 1 / 2 at very early times. On the other hand, Cornell et al 8 ' 9 ) have 
recently treated with the diffusion - limited reaction nA + mB^C 
for both homogeneous and inhomogeneous conditions under initially 
separated reactants. They have found that the global reaction rate 
decreases as t~ 1 / 2 at long - time regime, independent of n and m, and 
discussed that the upper critical dimension is d = 2 for the reaction - 
diffusion process. Very recently, Yen et al 10 ) also have studied for the 
asymptotic early - time scaling in the ternary A + 2B^C reaction - 
diffusion process with initially separated reactants. To our knowledge, 
it is of fundamental importance to measure extensively more compli- 
cated reaction such as multi - components reaction - diffusion process, 
even though the three - species reaction may occur rarely in the fields 
of physics, chemistry, and biology. 

In the present study, we present the decay process in a reaction- 
diffusion system with three - species. By a simple perturbation ex- 
pansion, we derive analytically the particle density and the global 
reaction rate before the crossover in the reaction - diffusion system of 
A + B + C — > D. The scaling exponent for the global reaction rate is 
also confirmed numerically before and after the crossover via Monte 
Carlo simulations. 

First of all, let us denote that pa(x, t), Pb(x, t), and pc(x, t) are the 
particle densities for three - species A, B, and C existing at a position 
x at time t. As we assume that three - species particles are distributed 
separately on the axis x, the special initial condition for these parti- 
cle densities is as follows : p^(x,0) = AqH(x), pb(x,0) = BqH(x), 
and pc(x,0) = Cq[1 — H{x)\. Aq, Bq, and Cq are the initial particle 
densities, and H(x) is the Heaviside step function. We assume that 
both Aq and Bq have a uniform mixing state on < x < oo. Since the 
mean field approximation may be applied to three - species system, 
A + B + C — ► D, the deterministic rate equation for pA(x,t) can be 
expressed in terms of 
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q^Pa(x, t) = D A V 2 p A (t) - Kp A (x, t)p B (x, t)p c (x, t), (1) 

where D A is the diffusion constant for one - species A. In the above 
equation K is the reaction rate which can be considered as small value 
in reaction - limited reaction. 

In particular, when we introduce the dimensionless particle den- 
sities a(x,t), f3(x,t), and j{x,t) from the relations of p A (x,t) = 
A a(x,t), pB(x,t) = B /3(x,t), and pc(x,t) = C j(x,t), eq.( 0) for 
a(x, t) can be rewritten as 

^-a(x, t) = d -^a(x, t) - -a(x, t)f3(x, t)~/(x, t), (2) 
at rj 

where the perturbation parameter do, r\, and u are denoted by do = 
(D A /D B D C ) 1/2 , V = 4), and u = K / {D A D B D C ) 1/2 . In the pertur- 
bation parameters Db and Dq are the diffusion constants for two - 
species B and C. For early - time regime, a perturbation theory can 
be developed to calculate on the decay process of three - species re- 
actants , provided that the reactive effect is small compared to the 
diffusion effect. For the sake of concreteness, if it is assumed to be 
extremely small (i.e. u<l) for the perturbation coefficient, the di- 
mensionless particle densities are expanded in series of powers of u> as 
a = En=i a n u n , P = En=i / W, and 7 = En=i lnOJ n - The solutions 
for a(£, r), r), and 7(£,t) are obtained that 

aK ' T) ^ ( ^r ) '' 5 «' T) ^ ( ^ > - (3) 



and 



7(£,t) = 1-«4(^) 1/2 ], ( 4 ) 



where the dimensionless time r and space variables £; are, respec- 
tively, given by r = tAoB Co(D A D B D c ) 1/2 , £ = x(A BoC ) 1/2 , and 



[x) = 4p fZ^ dxe x2 



In general, the global reaction rate R(t) also can be defined by 
R(t) = dxR(x, t), where the local reaction rate is the statisti- 
cal quantity given by R(x,t) = Ka(x,t)(5(x,t)'y(x,t). For simplicity, 
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we can get the global reaction rate from lowest order of perturbation 
terms. Hence, this is calculated as 

R(t) = A B Co{D a D b Dc) 1/2 uj{R +ujR 1 + ...) f 



in early - time regime. 

In order to confirm numerically the analytical result of the global 
reaction rate, we concentrate mainly on the Monte Carlo simulation 
in reaction - diffusion process of A + B + C — > D. It is supposed 
that three - species reactants are distributed randomly on two - di- 
mensional square lattice with the periodic boundary conditions. We 
also assume that the intermediate process existing concurrently the 
combined two - species reactants can be formed, when two - species 
reactants meet each other on the same lattice points. If two - species 
reactants meet the third reactant, these reactants react and leave im- 
mediately on lattice point. In this reaction - diffusion process, we only 
restrict ourselves to the case that the diffusion constant takes the same 
value for each reactant. 

In our reaction - diffusion system, the respective particle density 
of 30% for A, B, and C is distributed randomly on a square lattice 
having 200 x 200 lattice points. When our simulations are performed 
on 2 x 10 2 realizations, our reaction constant is 1/1500. As shown 
in Fig.l, we can directly observe the crossover for the global reaction 
rate from our simulation result. Therefore, it is found numerically 
from Fig. 2 that the scaling exponent for the slope before the crossover 
is 0.49 that is close nearly to 0.5 , the result of eq.( ||) derived from 
our three - species process. The scaling exponent for the slope after 
the crossover also is —0.54, as shown in Fig. 3. 

In conclusion, we have analytically derived the global reaction 
rate before the crossover in reaction - diffusion system of the form 
A + B + C — > D. It is really found from our simulation result that 
the global reaction rate is approximately propertional to t 1 / 2 at early 
- time regime, and to i~ x / 2 at long - time regime by scaling argu- 
ments. However, future work is in progress to extend to the system- 
atic method of renormalization field theory for multi - components 
reaction - diffusion process, and we also will attempt to investigate 
extensively on fractal lattice models. 11 ) 
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Figure Captions 



Fig.l. 

The global reaction rate R(t) verse time t for the reaction - diffu- 
sion process of the form A + B + C — > D. Our numerical simulation 
is done for 2 x 10 2 configurations on 200 x 200 square lattice with the 
periodic boundary conditions. 



Fig.2 

Plot of lnR(t) verse Int before the crossover. The slope of solid line is 
0.49. 



Fig.3 

Plot of lnR(t) verse Int after the crossover. The scaling exponent 
we obtained is approximately —0.54. 



6 



